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We present the studies of non- uniformly distributed vortices in Nb /Al mul- 
tilayers at applied field near parallel to film surface by using spin-polarized 
neutron reflectivity (SPNR) and DC magnetization measurements. We have 
observed peaks above the lower critical field, Hd, in the M- if curves from the 
multilayers. Previous works with a model calculation of minimizing Gibbs 
free energy have suggested that the peaks could be ascribed to vortex line 
transitions for spatial commensuration in a thin film superconductor. In or- 
der to directly determine the distribution of vortices, we performed SPNR 
measurements on the multilayer and found that the distribution and density 
of vortices are different at ascending and descending fields. At ascending 
2000 Oe which is just below the first peak in the M-H curve, SPNR shows 
that vortices are mostly localized near a middle line of the film meanwhile the 
vortices are distributed in broader region at the descending 2000 Oe. That 
is related to the observation of more vortices trapped at the descending field. 



As the applied field is sightly tilted (< 3.5°), we observe another peak at a 

smaller field. The peak position is consistent with the parallel lower critical 

field (Hci||). We discuss that the vortices run along the applied field below 

Hci|| and rotate parallel to the surface at Hdy. 
PACS numbers: 74.60.Ge, 61.12.Ha, 74.25.Ha 
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I. INTRODUCTION 



Vortices running parallel to surface in a thin superconducting film have been widely 
investigated theoretically and experimentally. As the field applied parallel to surface, the 
Bean-Livingston surface barriers significantly contribute to ingress and egress of vortices 
|T|J^. They also could do an important role in determining density and distribution of 
vortices, as the film thickness is comparable to the London penetration depth Unusual 
prominences in a M-H curve above Hd have been observed from superconducting films by 
using several different techniques, including electron tunneling 0, microwave absorption 
IP, resistivity 0, superconducting channel device SQUID magnetization torque 
magnetization [Q, and vibrating reed [|l^,|l^. Guimpel et al. |^ has suggested that the 
peak could be due to the vortex line transitions for spatial rearrangements with a model 
calculation of minimizing Gibbs free energy. The idea of minimizing free energy with vortices 
has been further developed by Brongersma et al. with Monte Carlo simulation . However 
the above techniques measure only the average magnetization or the result of vortex motions. 
It means that they can not determine the location of vortices with the mesurements alone. 



Furthermore experimental measurements from a thin YBa2Cu307_i. film |10| disgree with the 



model calculation of the free energy with vortices [|TT|. Therefore only a direct measurement 
on the location of vortices would clarify whether the peaks comes from the vortex line 
transitions. 

Since spin-polarized neutron reflectivity (SPNR) could detects spatial gradient of mag- 
netic field, it has been used to determine London penetration depth of high-Tc superconduc- 
tors as well as conventional superconductors [|l^ at a small field. In the regime of saturated 



field, SPNR has showed a capability to observe vortices |jT3|,|T^. A theoretical model calcu- 



lation by Han et al. [|13| has demonstrated that it could be able to measure the density and 
also the distribution of vortices. By using another advantage of SPNR that measures only 
parallel component of magnetization due to a scattering selection rule, Han et al. [|15| has 
demonstrated that vortices placed perpendicular to surface considerably contribute to the 
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magnetization measurements even at a small tilted field. 

In this paper, we introduce SPNR and magnetization measurements from the vortices 
in Nb/Al mulitilayers in Sec. II and present the results in Sec. III. A model calculation for 
minimizing Gibbs free energy is discussed in Sec. IV and the magnetization measurements 
under tilted fields are presented in Sec. V. In the last section, we summarize the main 
conculsions. 

II. EXPERIMENTAL DETAILS 

Nb(72A) /A1(20A), Nb(lOOA) /A1(20A) and Nb(130A) /A1(20A) mulilayers with repeating 
20 times respectively were deposited on Si substrates by direc-current sputtering under a 
base pressure of ~10~^mTorr and the Ar partial pressure of SmTorr. During the deposition, 
the power was apphed to a Nb target with 275 watts (297 voltages) and a Al target with 
200 watts (372 voltages) at which the deposition rate was 5.9 A/ sec for Nb film and 4.7 
A/sec for Al film meanwhile the substrate was placed at the ambient temperature. Tc of the 
Nb(72A)/Al(20A) multilayer, 7.25 ± 0.25 K, was determined by magnetization measurement 
at applied field 50 Oe. 

For magnetization measurements, a multilayer mounted on an extended sample holder 
was placed in a cryostat with the geometry where the film surface runs near perpendicular to 
a pick-up coil of the SQUID and zero-field cooled. Subsequently the field was applied. The 
tilt angle between the surface and applied field was controlled by shimming non-magnetic 
plastic pieces between the sample holder and sample. The tilt angle was reproducible within 
uncertainty, ±0.25". 

SPNR measurements from the Nb(72A) /A1(20A) multilayer were performed at POSYl 
reflectometer. Intense Pulsed Neutron Source, Argonne National Laboratory. The polar- 
ization efficiency ~97% and the instrumental resolution Aq/q = 0.053 were counted to 
analyzing the reflectivity data. The specimen was zero-field cooled and the field was applied 
parallel to the film surface at tilt angle < 0.5°. The polarized neutrons was reflected from 
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the multilayer with incident angle 0.45". 



III. DATA ANALYSIS AND RESULTS 

Magnetizations from the multilayers were measured for ascending field at 2 K, as shown in 
Fig. |I|. We observe the peaks (indicated by small arrows) from the multilayers. The peak po- 
sitions are dependent on the configuration of the films. The peak from the Nb(72A)/Al(20A) 
multilayer is more notable than the others. The detail magnetization studies will be dis- 
cussed in Sec. V. For more understanding of the interfaces of the Nb(72A)/Al(20A) multi- 
layer, specular x-ray reflectivity measurement was carried under atmosphere. Figure ^ shows 
the grazing angle x-ray reflectivity as a function of momentum transfer, q. A best fit (solid 
line) shows Nb(120±20A)/[Al(19±1.5A)/Nb(74.5±2.5A)]x20/Si and an extra layer at the 
air/Nb interface with thickness ~60 A and the x-ray scattering density half of Nb improves 
a best fit. A rms roughness at the air/Nb interface is ~15 A while the roughness of the Si 
substrate is ~2.2 A. Each Nb layer has the roughness ~6 A and each Al layer has ~2.5 A 
rough surface. 

SPNR measurements were performed on the Nb[72A]/Al[20A] multilayer to directly de- 
termine the density and distribution of vortices which run parallel to the surface. Figure |^ 
(a) shows neutron specular reflectivities measured as a function of q for spin up and spin 
down neutrons at 700 Oe and 2 K. The oscillation period is corresponded to the total film 
thickness ~2020 A. A best fit (solid line) shows that the top Nb layer thickness is 180 ± 
40 A and multilayer [Al(20 ± 2A)/Nb(72 ± SA)] x 20. However it could not reveal the 
buried interfaces as well as the x-ray reflectivity does because the neutron data were taken 
in a small- g region. The neutron reflectivity is basically consistent with the x-ray reflectivity 
within uncertainty and the configuration of the multilayer determined by neutron reflectivity 
measurement was used for analyzing the SPNR data. Figure |^ (b) shows reflectivity dif- 
ference between spin-up and spin-down devided by their average (AR/i?) that more clearly 
demonstrates magnetization contribution to SPNR. The solid line is a best fit without in- 
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eluding vortices. It showes that the London penetration depth (A^) of the multilayer is 1800 
± 200 A. At a high field, vortices will enter the superconductor and SPNR would see the 
vortices. 

Figure^ show AK/R under 2 K at ascending field 1500 Oe (a) which is just above Hd and 
2000 Oe (b) which is just below the first peak, as shown in Fig. |^. Below critical angle (g ~ 
0.013 A~^), the beam is totally refiected by the surface and the magnetization contribution 
to SPNR vanishes. For q > 0.018 A~^, the refiectivity differences among the models are 
relatively smaller meanwhile the uncertainty is growing up. Thus, data in g-range of 0.013 



- 0.0175 A^^ only were fit to a theoretical model ||13[. In the fit, only single parameter, 
density of vortices, was varying and the distribution is assumed to be mostly localized in 1 
row (dashed line), 2 rows (solid line) and uniformly distributed through the whole specimen 
(dotted line). 

The results of least x squares fits with counting the statistical uncertainty are summarized 
in table | for 1500 Oe. It shows that 1 row of vortex model is slightly better than the others 
and the vortex density is 30 ± 6 /im~^ which is corresponding to the distance of vortex- 
vortex ~1650 A. Table y shows the of fits for different models for 2000 Oe. The models of 
1 row and 2 rows are not distinguishable with basing on the least x^- However the magnetic 
field due to the vortices in a line interferes with the surface screening fields in AR./R. The 
interference effect is more clearly seen near the first peak of AK/R, q ~ 0.0145 in this 
system. It suggests that the model of 1 row is better than 2 rows. From the fit with 1 row 
(dashed line) the vortex density is found to be 45 ± 6 /im~^ that is corresponding to ~1100 
A of the average distance between two adjacent vortices. 

After cycling the field to 5400 Oe, SPNR measurement was also conducted from the 
multilayer at 2000 Oe, as shown in Fig. ^. The lines present a best fit with assuming 
different distributions of vortices. A goodness fit of solid line strongly suggests that the 
vortices stay in 2 rows, 1/3 and 2/3 of the total film thickness and the density of vortices 
56 ± 2 yum^^ The results of fits are summarized in table The fits show that 2 rows of 
vortex model is about twice better fit than 1 row and uniform distribution. From the best 
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fit with 2 rows, the nearest vortex distance is calculated to be ~1110 A with assuming a 
triangular vortex lattice, (b) shows the neutron scattering density profile of magnetic as 
well as nuclear potential that is corresponding to the fit with 2 rows in (a). The inset shows 
a vertical expansion. 



The best fits of the SPNR data are summarized in table The analyses show that 
vortices stay in a single row at ascending fields, 1500 and 2000 Oe, however they would like 
to spread into 2 rows at the cycled field 2000 Oe. More broadening at the cycled field could 
be explained in terms of the trapped vortices by the surfaces. The parallel magnetization 
(M||) is calculated by using the density and distribution of vortices determined by SPNR 
with a theoretical model |]T5[. The M|| shows that is unable to distinguish the differences 
among the fields due to large uncertainty. However they are certainly negative values even 
at the cycled field. We could directly compare those to the SQUID data. As shown in the 
table, the magnetizations of M|| for ascending field are comparable to M directly converted 
from SQUID data within uncertainty however there is a big difference at the cycled field. It 
suggests that M of the cycled field is contributed by not only M|| but also vortices running 
perpendicular to the surface. That mechanism in which vortices enter the superconducting 
film parallel to the surface for ascending field and rotated out of plane during reducing the 



field has been observed in a Nb superconducting film [15 



Although we could obtain reasonable results with the model, 1 row or 2 rows, of the 
vortex distributions, a Gaussian distribution could be an alternative choice for the vortices 
in a superconducting film. Figure ^ shows the best fits with a Gaussian distribution of 
vortices and table |VI| summarizes the results of the fits. The model of Gaussian distribution 
absolutely improves the fits, particularly 1500 and ascending 2000 Oe. At 1500 and ascending 
2000 Oe, the full width at half maximum (FWHM) is basically the same however it shows 
broader at the cycled 2000 Oe. The density and M|| are very comparable to the model of 1 
row or 2 rows. Because of a limit of the SPNR sensitivity to this system, it is not able to 
distinguish between a Gaussian distribution and 1 row or 2 rows. However it is clear that 
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the vortices have a same distribution at ascending fields and they are differently distributed 
at the cycled field. 



IV. MODEL CALCULATION AND DISCUSSION 

For more theoretical understanding of vortex distribution in a superconducting film, 
we have developed a model calculation for minimizing Gibbs free energy. As N number of 
vortices enter in a superconducting film, the total free energy of the superconductor could be 
calculated by using a London approximation p| with simply assuming that it is an isotropic 
superconducting film and impurity in the film is negligible. 

G = Go + ^f dv{$-i2BL + Bv)} (1) 

Zfig JVS ^ 

where Gq is free energy for the system without vortices, /Iq is permeability in vacuum, Vs is 
volume of superconductor, $ is vorticities, Bl is surface screening field and By is magnetic 
field due to vortices including their images. The vorticities are defined to be, 

k=N 

$ = $0 ^ 5(f -ffc)£ (2) 

k=0 

where $0 is a vortex fiux quantum, 2.067 x 10^ GA^, is location of k*'^ vortex and vortices 
are oriented in x-axis. As assumed that the applied field is exponentially decayed from the 
surface, the surface screening field is 
^ f cosh(z/ Xr) 1 „ 

where t is film thickness and applied field is along x-axis. If the vortices are localized in lines, 
•r k = Zk z + kl y in Eq. (Q) where I is average distance of adjacent vortices in ^-direction and 
the the spatial magnetic field due to the vortices will be 

^ k=N/2 „ r J(z -nt~ i-lYzkf + {y- km ^ 



fc=_7v/2 n=-oo 

where Kq is a modified Bessel function of the first order. For simplifying the free energy 
calculation, previous studies PJl8[] have used an approximation which is valid only for a 
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limit of ^ a film thickness. Since the approximation is not applicable for our system (A^ 
< t), we have to count term by term in the summation of Eq. (^). 

For the free energy calculation, we need two characteristic lengths, London penetration 
depth and coherence length. A^ ^ 1800 A for the Nb(72A)/Al(20A) multilayer was deter- 
miend by SPNR. The coherence lenght could be obtained by measuring Hd. As the field 
applied parallel to the surface, the low critical field can be estimated by using the London 
theory With assuming that a vortex first enters a thin superconducting film as the 

free energy is zero at z where the vortex is placed, the lower critical field is 

$n 1 



Hci 



X 



27rAi 1 - cosh{z/\L)/cosh{t/2\L) 



We assum that vortices first enter the superconductor, as the free energy is zero at the 
surface, e.q., z = t/ 2- ^in Eq. (||) [0. Based on the magnetization measurement, Hd = 
1200 ± 200 Oe and ^ is calculated to be ~113 A. Although the Nb(72A)/Al(20)A muhilayer 
is an anisotropic superconductor, an anisotropy of a Nb(100A)/Cu(100A) multilayer, iy/iz, 
was 1.23 0], where and are the coherence lengths of in-plane and out of plane respec- 
tively. As comparing configuration of our specimen to Brongermas's one, the asumption of 
an isotropic superconductor for the calculation would not be seriously wrong. 

Fig. 1^ shows the minimum free energy calculation. For the calculation, it is assumed 
that vortices are localized in a central line {t / 2) (dashed line) and two lines {t / 3 and 
2t I 3) (solid line), as shown in the inset. At a given field, the minimum free energy was 
determined by varying the density of vortices only. At small fields, the free energy of the 
system for 1 row is smaller than for 2 rows whereas above 2200 Oe it is reverse. It means 
that the vortices more likely stay in 2 rows than 1 row for H > 2200 Oe. This calculation 
agrees well with the first peak in the magnetization measurement, as shown in Fig. |l] top 
and strongly suggests that the peak is due to a vortex line transition from 1 row to 2 rows. 
The inset at upper right corner shows magnetization which corresponds to the minimum 
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free energy. There is the second transition (2 3) at ~4000 Oe that is also consistent with 
the second peak in the M-H curve. 

The calculation shows that vortex density of 40 fim~^ at 1500 Oe and 53 /im~^ at 2000 
Oe will satisfy the condition of minimizing the free energy with 1 row of vortices in this 
system. From SPNR measurements, it is found to be ~30 /im~^ at 1500 Oe, ~45 /im"^ 
at 2000 Oe and ~56 yum~^ at cycled 2000 Oe. The smaller vortex densities at ascending 
fields can be understood in terms of the surface barriers because the minimizing free energy 
calculation does not count the barriers. Since the vortex densities at ascending fields are 
smaller than the maximum vortex density for 1 row, the vortices could stay in a central line 
however at the cycled field, the vortex density found by SPNR is slightly higher than the 
the maximum vortex density for 1 row. Thus they can not stay in a single row. It means 
that the vortex line transition fields can be shifted to a lower field for descending field, as 
the contribution of the surface barrier is important. It has been experimentally observed by 
J. Sutton et al. 0. 

V. MAGNETIZATION AT SMALL TILTED FIELD 

As the field is tilted with a small angle, the delay of vortex entrance due to the sur- 
face barrier vanishes because the perpendicular component of applied field helps the vortex 
ingress at a small field. Figure |^ shows the magnetization at different tilt angles. There is 
a peak at H ~ 2250 Oe from Nb(72A)/Al(20A) film without depending on the tilt angle 
{9 < 3.5°). However we find another peak at ~950 Oe of the tilted field. It is the first 
observation that the first peak is missing at no-tilted field. One could easily overlook that 
the first peak does not come from the vortex line transition. Since it is a thin film, vortices 
can not enter the superconductor until the applied field is stronger than the surface barrier 
field. At a tilted field, however, vortices can overcome the surface barrier under support 
by the perpendicular component of the field even with a small field. Therefore the vortex 
line transitions could be missing at no-tilted field. This scenario could work only for that 
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the surface barrier field is larger than the vortex line transition fields. That is not this case 
because the measured surface barrier field for the specimen is about 1200 Oe and SPNR 
measurement shows 1 row of vortices at ascending fields of 1500 and 2000 Oe. Also the 
calculation for minimizing the free energy suggests that the first vortex line transition may 
occur at ~2200 Oe. Thus, we suspect that the first peak at the tilted field comes a vortex 
line transition. 

Using Eq. (|^), Hd of this system under the equilibrium condition where vortices first 
enter the superconducting film as the free energy is zero at the middle of the film, e.q., z = 
in Eq. is calculated to be ~850 Oe. That is very comparable to the first peak at the 
tilted field. It might suggest that the peak is due to the lower critical field which is paralle 
to the surface (Hcij|). At a tilted field {6 > 1° for this system), the first vortex entrance is 
determined by Hci± instead of Hci|| because for a thin film superconductor, Hci_l is much 
smaller than Hciy [0. For B.sin9 (B.±) > iiciL, vortices could enter the superconductor and 



might stay along the applied field due to the dragging force from the applied field. However 
the free energy still does not allow the vortices running parallel to the surfaces. When the 
applied field (HcosO) is increased with passing beyond Hci|j, vortices sightly rotate to parallel 
to the surfaces, connect the pieces, and make long threads. In this scenario, however, the 
magnetization might be just a little changed because the field is applied with a small tilt 
angle. However we find that a large demagnetization effect enhances the change by more 
than two orders of magnitude because of the thin film geometry. The demagnetization of 
the multilayers is discussed in detail below. 

Figure ^ (b) and (c) show the magnetization from the multilayers at 4.5 K. The magne- 
tization from Nb(100A)/Al(20A) also shows that the first peak is missing at no-tilted field 
whereas it is observed at 600 Oe of a tilted field. That agrees well with Hdn = 600 Oe 
calculated by using Eq. (ID (z = 0), and = 1800 A and ^ = 113 A. The second and third 
observed peaks, ~1400 and ~2450 Oe, can also be compared to the vortex line transitions, 
1^2 (1450 Oe) and 2— > 3 (2650 Oe) respectively estimated by the model calculation. That 
is another evidence that the first peak at a tilted field comes from the vortex rotation instead 
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of the vortex line transition. At 4.5 K, two peaks are observed at 845 and 2125 Oe from 
Nb(72A)/Al(20A). The first peak near 845 Oe appears without depending on the tih angle. 
It implies that the superconductor is soft and the surface barrier is negligible at the temper- 
ature. Since the peak positions are not very sensitive to the temperature that is consistent 
with previous measurements , we might also conclude that the first peak comes from the 
vortex rotation at B.ci\\- The peaks from the Nb(72A)/Al(20A) and Nb(100A)/Al(20A) mul- 
tilayers were carefully measured at different temperatures and tilt angles and summarized 
in table Table shows that the vortex line transitions (2nd and 3rd peaks) occur at 
smaller fields for a thicker film superconductor because vortices have relatively more space 
to the surfaces in a thicker film and rearrange for spatial commensuration at a smaller field 

iiHl- 

Demagnetization factors of the multilayers are quantitatively determined in the Meissner 
regime. In the regime, slopes of the magnetization were measured as a function of tilt angle, 
as shown in Fig. ^ The solid lines are a best fit to data with a model for the demagnetization 
of a thin superconducting film |T5|. The demagnetization factors of Nb(72A)/Al(20A) are 
found to be 0.9986 ± 0.0011 at 2 K and 0.9935 ± 0.0007 at 4.5 K while it is 0.992 ± 0.0034 
for Nb(100A)/Al(20A) at 4.5 K. Those can be compared to theoretical calculations, 0.993, 
0.9942 and 0.9934 respectively for our sample geometries |2^. For H < B.ci\\, the perpendic- 
ular magnetic field due to vortices is n^effsin6/{l — N) where n is number of vortices per 
unit area, $oe// is the effective flux due to the images and N is the demagnetization factor. 
As the vortices rotate parallel to the surface at Hdn, the perpendicular magnetic field will 
vanish. As assumed that a single vortex per /im^ rotates parallel to the surfaces at the tilt 
angle 2.5", the difference of a few Gauss for our sample geometry contributes to the SQUID 
magnetization measurement. 
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VI. CONCLUSIONS 



Vortex pinning is an important subject for practical application of current transportation. 
With only thin film geometry where vortices place parallel to the surfaces, the vortex-surface 
interaction could be studied. We performed SPNR measurements on a Nb/Al multilayer to 
directly determine the density and distribution of vortices. SPNR shows that the vortices in 
a Nb/Al multilayer are localized near the central line below the first peak in the M-H curve 
at the ascending field meanwhile vortices have more broadening of the distribution at the 
descending field. The SPNR measurement is consistent with the magnetization measurement 
and also a model calculation of minimizing Gibbs free energy. As the field is apphed with an 
angle, another peak which is missing at no-tilted field is first time observed. Comparing the 
peak to a model calculation of Hciy, we conclude that the first peak at a tilted field comes 
from the vortex rotation instead of vortex line transition at Hciy. 
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FIG. 1. shows magnetization measured from Nb/Al multilayers at 2 K. The data of Nb(100 
A)/A1(20 A) and Nb(130 A)/A1(20 A) were shifted down along the y axis for clarity. The smah 
arrows indicate the peaks and the big arrows points the direction of the measurements. 

FIG. 2. X-ray reflectivity was measured from the Nb(72A)/Al(20A) multilayer at atmosphere. 
Dotted line is data and solid line is a best fit. 

FIG. 3. (a) shows grazing angle reflectivities for spin up and spin down neutrons from a 
Nb(72A)/Al(20A) multilayer measured as a function of q at 700 Oe and 2 K. The solid line is 
a best fit. (b) shows AR/R (described in the text) obtained from the data in (a). The solid line 
is a best fit and shows that the London penetration depth of the multilayer is ~ 1800 A. 

FIG. 4. shows AR/!R measured from a Nb(72A)/Al(20A) multilayer at ascending 1500 Oe (a) 
and 2000 Oe (b) under 2 K. The lines are a best fit with assuming different distributions of vortices. 
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FIG. 5. (a) shows AR/R measured from a Nb(72A)/Al(20A) multilayer at cycled field 2000 Oe 
under 2 K. The lines are a best fit with assuming different distributions of vortices, (b) shows the 
neutron scattering density profile that is corresponding to 2 rows of vortices in (a). The inset is a 
vertical expansion. 

FIG. 6. shows AR/R at 1500, 2000 and (cycled) 2000 Oe and the solid lines are a best fit with 
a Gaussian distribution of vortices. 

FIG. 7. shows minimum Gibbs free energy as a function of applied field with assuming 1 row 
(dashed line) and 2 rows (solid line) of vortices. The assumed locations of vortices are shown in 
the inset at lower left. The inset at upper right shows the magnetization corresponding to the 
minimum free energy. 

FIG. 8. Magnetic moments were measured from Nb/Al multilayers at different tilt angles and 
different temperatures. The arrows in (a) indicate the direction of measurement. Data in (b) and 
(c) were taken along the same direction of measurement in (a) . The lines are a guide to the eye. 

FIG. 9. shows the slopes of magnetic moment below Hd as a function of tilt angle. The solid 
lines are a best fit with a model [15]. The data (open circle and solid triangle) are shifted down 
along the y axis for clarity. 



TABLE I. Least-x^ fit at 1500 Oe 





1 row 


2 rows 


Uniform 




2.788 (30 Aim^2) 


2.974 (30 /Lim-2) 


3.499 (35 /Lim-2) 
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TABLE II. Least-x^ fit at 2000 Oe 





1 row 


2 rows 


Uniform 


,2. 
X ■ 


z.y4o (40 /iiii j 


z.yix) (40 /im j 


o.Uo4 (dU /im j 




TABLE III. 


Least-x^ fit at cycled 2000 Oe 






1 row 


2 rows 


Uniform 




1.855 (50 


0.969 (56 


1.498 (72 /im-2) 




TABLE IV. 


Results of SPNR measurements 




H (Oe) 


Distribution 


Density per fim? My (G) 


M (G)^ 


1500 


1 row 


30 ± 6 -51 ± 17 


-66.6 ± 7 



2000 1 row 45 ± 6 -57 ± 17 -62.8 ± 6 

2000 (cycled) 2 rows 56 ± 2 -42.5 ± 5 83 ± 8 
'^Magnetization directly converted from SQUID data in Fig. |8| (a) open circle with no adjustable 
parameter. 

TABLE V. Least-x^ fit with a Gaussian distribution of vortices 

H (Oe) FWHM (A) vortex density (fim-^) My (G) 

1500 2.056 510 ±51 33 ± 3 -47 ± 9 

2000 2.077 530 ± 29 47 ± 3 -58 ± 9 



2000 (cycled) 0.774 


700 ± 23 


55 ± 2 


-42 ± 6 


TABLE VI. Peak positions(Oe) 


Configuration of multilayer 


1st 


2nd 


3rd 


Nb/[Al(20A)/Nb(72A)]x20 
Nb/[Al(20A)/Nb(100A)]x20 


950(2K, 4.5K) 
650(4.5K) 


2250(2K, 4.5K) 
1400(2K, 4.5K) 


3900(2K) 
2450(2K, 4.5K) 
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m (emu) 



1 


1 


1 


1 


1 
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o 
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o 
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o 




o 


oo 


cn 


4^ 




o 




Reflectivity 




2 - 




0.012 0.016 0.020 
q (A-l) 



Uniform distribution 




0.10 



0.05 



Uniform distribution 

1 row vortices 

2 row vortices ^ 







-0.05 - 



-0.10 




T = 2 K 
H = 2000 Oe 



- (a) 



0.012 0.014 0.016 0.018 



1.0- 



(b) 



0.5- 








500 



1000 
z (A) 



1500 2000 




0.010 0.015 0.020 

q (A-') 



(a.u.) 

I I 



cn 
o 



ho 
cn 



o 
o 



o 

cn 




0.010 



(a) 



Nb(72A)/AI(20A) 



0.005 




CD 







-0.005 




T = 2 K 



tilt angle 

^2.5° 

-e-0.0° 



500 



1500 
H (Oe) 



2500 



-0.0005 



-0.0015 



-0.0025 



Nb(100A)/AI(20A) 
T = 4.5 K 




*A m ^''^ angle 
•^•^ ^1.6° 
-A- 0.9° 
-•-0.2° 




tilt angle 



^2.75^ 
-A- 1.4° 
-•-0.0° 

T = 4.5 K 
Nb(72A)/AI(20A) 



500 10001500 500 10001500 
H (Oe) H (Oe) 



